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Abstract 

Consider a small sample coupled to a finite number of leads, and assume that the total 
(continuous) system is at thermal equilibrium in the remote past. We construct a non- 
equilibrium steady state (NESS) by adiabatically turning on an electrical bias between the 
leads. The main mathematical challenge is to show that certain adiabatic wave operators 
exist, and to identify their strong limit when the adiabatic parameter tends to zero. Our 
NESS is different from, though closely related with the NESS provided by the Jaksic-Pillet- 
Ruelle approach. Thus we partly settle a question asked by Caroli et al in 1971 regarding the 
(non)equivalence between the partitioned and partition-free approaches. 
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1 Introduction 



1.1 Generalities 

This paper deals with the rigorous construction of adiabatic non-equilibrium steady states for 
mesoscopic systems which initially are fully coupled (or 'partition free') and at thermal equilibrium 
[9"lll4j. The initial equilibrium state is broken down by slowly turning on an electrical bias between 
leads (i.e. inserting a d.c. battery), which in a certain way can be seen as slowly changing the 
chemical potentials of the leads coupled with the small sample. 

In contrast with the above described partition-free setting, the 'partitioned procedure' is the one 
in which one starts with several decoupled reservoirs, each of them being at different equilibrium 
states. Let us assume for simplicity that they are in grand canonical Gibbs states having the same 
temperature but different chemical potentials. Then at t = they are suddenly joined together 
with a sample, and the newly composed system is allowed to freely evolve until it reaches a steady 
state at t = oo. From a mathematical point of view this approach is by now very well understood, 
see for example [TJ [TH1 [33J HS1 IS] and references therein. One can allow the carriers to interact in 
the sample [TS], and the theory still works. Note that even if we choose to turn on the coupling 
between the reservoirs in a time dependent way, the result will be the same |15j . 

One can ask which approach is more physical; here is a quote from a paper by Caroli et al 
[8] from 1971 -maybe the first very influential paper on the subject- who came with the following 
observation about the partitioned procedure: One might raise a major objection to the above 
procedure; it amounts to establishing first the dc bias, and only later the coupling between the 
barrier and the electrode. Physically, it is the reverse that is true; the transfer matrix elements 
are always there, and the dc bias is established afterwards; it is not obvious that the corresponding 
limits can be interchanged. 

The major achievement of our current paper is that we can now construct an adiabatic NESS in 
the partition free setting; let us explain how. The leads are already coupled with the sample, and 
at t = — oo the full system is in a Gibbs equilibrium state at a given temperature and chemical 
potential. Then we adiabatically turn on a potential bias Vx(vt) between the leads, modeling 
in this way a gradual appearance of a difference in the chemical potentials (here x(~°°) = 0; 
x(0) = 1 and -q > is the adiabatic parameter). The final bias V does not need to be small; our 
results are beyond the linear response theory. The statistical density matrix p^it) is found as the 
solution of a quantum Liouville equation, with the initial condition at t = — oo given by the global 
Gibbs state. 

In Theorem 11.51 we show the existence and compute the strong limit p ac i := lim^o Pr)(f). The 
limit is t independent, and contains - as in the partitioned procedure- two contributions: one 
from the discrete, and one from the continuous subspaces. Note that we do not have to take the 
Cesaro limit in order to insure convergence for the discrete part. The adiabatic limit takes care of 
the oscillations. The price we pay is that we need to demand that the point spectrum of certain 
Hamiltonians only consists from finitely many discrete eigenvalues. Most probably this condition 
is too strong, and getting rid of it remains an interesting open problem. 

Even though the stationary density matrix of the partitioned procedure has a similar structure, 
it is different from the one we construct here. A careful comparison will be given elsewhere. 

A future problem is to investigate the charge current and establish Landauer-Biittikcr type 
formulas [6j [7j [2 |3j [14] in the partition free setting with a continuous model, and without the 
linear response approximation. In fact this was the starting point of a number of remarkable 
physical papers, see for example [17], [23], [4]. A first mathematically sound derivation of the L-B 
formula on a discrete model and under the linear response approximation was obtained in |12| and 
further investigated in [T3] . In QT] we significantly improved the method of proof of [12] , which 
also allowed us to extend the results to the continuous case. 

Another challenging open problem is to extend the formalism in order to accommodate tran- 
sient regimes (see (25] [26] [27] [14] and references therein), and locally interacting fermions [34l 135] . 

Finally, we want to stress that some of the technical conditions which we impose for our 
model (like smoothness of boundaries and potentials, working with only two parallel leads) can 



2 



be relaxed. We chose though to work under stronger conditions in order to give shorter proofs for 
certain spectral and asymptotic completeness results, thus making the paper rather self-consistent. 
In this way, the number of generic assumptions is kept to a minimum. 



1.2 The model 

Take two identical semi-inhnitc cylinders and couple them smoothly through a finite domain. The 
cylinders will model the leads, while the connecting domain will represent the region where the 
interesting physics takes place. The total configuration space £ is a subset of R d+1 with d > 0. In 
order to simplify presentation, we will assume that £ is cylinder- like, which means that for each 
value of the longitudinal coordinate x\\ £ M. the transverse coordinate belongs to a bounded 

*\\) 



cross-section 2)(xn) C M. d . Again for the sake of simplicity, we assume that the boundary 



£:=£>£ (1.1) 

defines a regular C°°-surface embedded in 

Let us start with the description of the configuration space associated to one of our d+1 
dimensional leads, namely the left one. Let a > 0. We let X := (— oo, —5) model its longitudinal 
dimension. Then we assume that: 

£n{X_ x R d } =:X_ x S, 

where the transverse section D C M d is supposed to be a bounded and simply connected open set 
with a regular C°°-boundary <92). Thus the configuration space of the left cylinder is modeled in 
a natural way by the set X_ x 2). Similarly, if X + := (a, oo), the configuration space of the right 
cylinder is modeled by X + x D. 
Now define: 

C:=£n{[-a,a] xK d }. (1.2) 

Thus the small sample is contained by a bounded and simply connected set C C W l+1 which is 
smoothly glued to the two leads. With these notations, the one particle configuration space can 
be decomposed as: 

£ = (X_ x ©) UCU (X+ x S). (1.3) 

When we refer to the " coupled system" , we mean that there are no internal walls between the 
sample and leads. A particle will be free to flow inside the system, and to pass from one lead to 
another via the sample. But it is not allowed to get out of £. 

Now let us introduce the one particle Hamiltonian of the coupled system. In the sample C we 
assume the existence of a potential iv € Cq°(C), which will be considered positive without loss of 
generality. The kinetic energy of a particle living in L will be modeled by the Laplace operator 
—Ad with Dirichlet boundary conditions on dC and having the domain Hd(L) := Hq (L) n H 2 (L) . 
Thus the one-particle Hamiltonian is of the form: 

H:=-A D +w, (1.4) 

with the same domain. 

Regarding the spectral properties of H , we will prove in Lemma l3.14l that its singular continuous 
spectrum is absent. We will assume that the pure point spectrum consists of discrete and finitely 
many eigenvalues: 

<7p P {H) = a disc (H), #a pp (H) < 00. (1.5) 

Remark 1.1. This assumption means in particular that we do not allow embedded eigenvalues in 
the continuous spectrum. To our knowledge, sufficient conditions to guarantee this property are 
not known in dimension d+1 > 2. 
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Let H := L 2 (L), and let a > a. Define 

L_ :=LD{(-oo,-o) x T>}, L+ := L n {(a, oo) x £>}, C := C H {(-a, a) x M d }. (1.6) 
We introduce three orthogonal projections: 

II : H -> H_ := L 2 (L_), n + : ft -> ft + := L 2 (L+), 

n :H^H :=i 2 (C). (1.7) 

Note that C is completely included in the open set C, and L± are "shorter" than the corresponding 
leads. 

1.3 The state and the Liouville equation 

Wc only work at the level of density matrices. In the remote past t — > — oo the electron gas is 
at equilibrium at a temperature T > and a chemical potential [i, moving in all the volume L. 
The gas is described by a quasi-free state, having as two-point function the usual Fermi-Dirac 
equilibrium density matrix operator: 

:= 1 + ei L, )/kT - (1-8) 

The system is driven out of equilibrium by slowly turning on an electric bias 

V = V-U- +V+U+, (1.9) 

where v± are real constants. We want to introduce the bias adiabatically with an adiabatic 
parameter 77 > 0, as a time-dependent potential V^t) := xivtyV- O ne should have in mind 
x(t) = e , but only a few abstract properties of this function are really needed, namely: 

0<x(*)<l and X '(t)>0 if t < 0; X (0) = 1; (1.10) 

x^'lei 1 ^-). 

We will also need to consider the 'bias' with a fixed coupling constant k £ [0, 1]. We introduce a 
family of operators: H(k) 

K{k):=H + kV. (1.11) 

Epp(A) and E ac (A) will denote respectively the projector on the pure point and absolutely continu- 
ous spectral subspace of the self-adjoint operator A. In Lemma l3.14l we will prove that the singular 
continuous spectrum of K (ac) is empty. We now make the following assumptions concerning the 
point spectrum: 

Hypothesis 1.2. 

1. Vk G [0, 1] the Hamiltonian K(k) has no eigenvalues embedded in the continuous spectrum; 

2. dimE pp (K(n)) = N < 00, G [0, 1], (t pp (K(k)) = {e 3 (n)}f =1 ; 

3. min {dist (a^fK (k)), a ac (K (k)))\ > d > 0. 
«e[o,i] 

In order to simplify the presentation, we will only work with N = 2 and adopt an extra 
assumption: 

Hypothesis 1.3. 

The eigenvalues {£j{n)}je{i,2} (which are real analytic functions of k 6 [0,1]) can cross at 
most at one point kq £ (0, 1). This ko corresponds to some unique to < where x(*o) = K o and 

x'(h) > 0. 
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The time dependent Hamiltonian will be 



K( X (vt)) := H + xivW (1.12) 

having the constant domain equal to the domain of i.e. IHd(L). The evolution defined by the 
time-dependent Hamiltonian K{x{ r l^)) is described by a unitary propagator W v (i), solution of the 
following Cauchy problem: 

" id t W v (t) = K( X (vt))Wr,(t) 

W v (0) = i, [LA6) 

for teR. For any ij > 0, the family {K(x(rit))}t^m consists of self-adjoint operators in H having a 
common domain equal to Hd(L) and strongly diffcrcntiable with respect to t e M with a bounded 
self-adjoint norm derivative d t K{x{ , r]t)) = rf^{rft)V. 

Now using well known results quoted in [311 Th. X.70] we easily obtain that the problem 
(I1.13[) has a unique solution which is unitary and leaves the domain IId(L) invariant for any 
t £ R. Moreover, its adjoint satisfies the equation: 

id t w;(t) = ~w;(t)K{x n {t)). (i.i4) 

The object we are interested in is the time evolved density matrix p v (t) which must be a 
solution of the Liouvillc equation, starting from the initial value p eq (H) at t — > — oo: 

id t p v (t) = [K(x{r]t)),p v (t)], n - limpet) = p eq (H). (1.15) 

t— > — oo 

In the remaining part of our paper we will show that the unique solution p^t) of the Liouville 
equation has a strong limit when 7/ \ 0, and compute it. In particular, we will see that the 
adiabatic limit is t independent. 



1.4 The main result 

In order to formulate our main result we need to define some new objects. First, we introduce 
the decoupled Hamiltonian obtained from H by introducing Dirichlet walls where the bias is 
discontinuous = ±a). Remember that the decomposition (|1.6[) depends on a, and the walls 

o 

are inside the leads. Let Ac be the self-adjoint Laplace operator defined in L 2 (L) with Dirichlet 

o o o o 

conditions on dL± U dC; we have = A^.- ffi ^-d.o ® ^d,+ , where their domains are denoted 
as follows: 

Bb(L±) := H^(L ± ) n H 2 (L ± ), U D (C) := H%(C) H 2 (C), 
U D {L) := Hd(L_) M D {C) © M D (L+). (1.16) 

Let us note that due to the cylindrical symmetry of the regions L± where the bias is piece wise 
constant, we can write 

A D ,± = i± 8> 1 + 1 ® £v (1.17) 

with £x> the Laplacean on the bounded domain T> cM. d with Dirichlet conditions on the boundary 
dT> 1 and l± the operator of second derivative on X± with Dirichlet condition at ±a. The decoupled 
one particle Hamiltonian will be: 

H:=-A D +w, (1.18) 

o 

which is self-adjoint on the domain Hd(L), having Dirichlet conditions on dL± U dC. As in the 
coupled case, we need to consider the bias with a fixed coupling constant k g [0, 1] and define 

o o 

K{k) := H + kV. In order to formulate our main theorem we need the following lemma: 

o 

Lemma 1.4.The stationary wave operator So associated to the pair {K(l), K(l)}: 

S := 8 - lime isK ^e- is ^E ac (K(l)), 

s\ — OG 
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exists and is a unitary operator from E ac (K(l)yH to £' ac (A'(l))H. Moreover, the singular contin- 
uous spectrum of K(k) is empty for all k £ [0, 1]. 

And here is the main result: 
Theorem 1.5. The adiabatic limit of the density matrix exists in the strong operator topology on 
B('H), is independent oft and given by: 

N 

Pad := s - lim p v (t) = E oPeq (H)E* + V p eg (£ 3 (0))^(tf(l)), (1.19) 
"NO ^ 

where {£j(0)}jLi are the eigenvalues of H = K(Q) in ascending order, while {E {K{Vf)}^ =l are 
the eigenprojections of H + V = K(l) obtained by analytically continuing {Ej{K{K)y\^ =l from 
k = to k = 1. 

Remark 1.6.. Even though Lemma lL4l is not surprising, its proof is not straightforward. 

We also note that the adiabatic limit p a d commutes with K(l) = H + V, but it is not a 
function of K (1). Even though Pr)(t) is a solution of a Liouvillc equation involving operators with 
no internal Dirichlet boundaries at ±a, the limit p ac i is expressed with the help of a comparison 

o 

operator H + V, depending on a, and which appears naturally in the proof. 

We will assume N = 2, but the result holds true for any finite N. An interesting open problem 
is to study the case N = oo and when the eigenvalues can enter the continuous spectrum while 
k grows from to 1. Another interesting situation is the one in which we have a degeneracy at 
k = 0; this situation is related to the Gell-Mann and Low theorem for degenerate unperturbed 
states [30] . 

1.5 A useful expression of the density matrix 

Before actually starting the study of the adiabatic limit, let us very quickly show that (|1.15j) has 
a solution, which can be put into a form which is particularly convenient for taking the adiabatic 
limit. 

Define the unitary adiabatic wave operators 



lu v :=n- lim W* (i)e~ ltjH , w* = n- lim e ltn W v (t). (1.20) 

t— y — oo t— > — oo 

They converge in norm due to the following estimate (s < t): 

\\W;{t)e- itH - W;{s)e- isH \\ < J* ±{ W ;{r)e-^ H } dr < \\V\\ j\(vr), (1.21) 

where we use that \ £ ^(R-). Then by direct computation we can prove that the operator 

Pv (t) := W v (t)u vPeq (H)u;w;(t) (1.22) 
solves the Liouvillc equation. It also obeys the initial condition because we can write: 
= lim - e-* tH {e ltH W v (t)}Lo vPeq (H)u:; {W*(t)e"^} e itH \\ 



= lim \\ Pn {t)~e- ttH p eq {H)e itH \\ = lim \\p v {t) - p eq {H)\\ . (1.23) 

The above solution can be rewritten as: 

p v (t)=W v (t) Pn (Q)W;(t), (1.24) 

where 

^(0) = Ur, Peq (H)u*. (1.25) 
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Now let us show that it is enough to prove (|1.19l) for t = 0. Indeed, once this formula is proved 

for t = it shows that the strong limit of p v (0) when r} \ is commuting with K(l) = H + V. 

It is elementary to check that W v (t) and W*(t) converge in norm to e ~ ltK ^ and respectively 

e itK(i) w } ien r] \ (with t fixed). Since e ±ltK ^ commutes with s — lim p„(0) it follows that the 

v\o 

adiabatic strong limit of p v (i) must also exist and equal the r.h.s of (|1.19|) . 

Moreover, due to the fact that the limits in (|1.20j) are in operator norm, it is easy to show that 
we have the identity: 

Pr,(0) =n- lim W*J S )p eq {H)W v (s). (1.26) 

s— > — oo 

It is important to note that the above norm limit is not uniform in rj, and this is the reason 
why the adiabatic limit is not straightforward. Formula (jl.26p will be the starting point in what 
follows, and we will be interested in computing the double limit: 




(1.27) 



2 A road map of the proof of the adiabatic limit 

Since our proof of the adiabatic limit is quite long, in this section we will give a list of technical 
results leading to it and postpone their proofs for the next sections. 

The two terms of (|1.19[) are coming from different spectral subspaces of H + V: the first one 
from the absolutely continuous spectrum, and the second one from the discrete spectrum. 

In Lemma 13.141 we will prove the absence of singular continuous spectrum for K(n), thus we 
can consider the orthogonal decompositions 

H = E ac ( K )H® { © Eji^uX, (2.1) 

where E ac {n) := E ac (K(n)) and Ej(n) := Ej(K{n)). Let us remark here the important fact that 
due to the Rellich Theorem (Theorem 11.61 in [20]) we can choose the eigenprojections of K(k) 
to be real analytic functions of n on the interval [0,1]. Then we can write 

w;( s ) Peq (H)w v (s) - w;( s ) Peq (H)E ac (o)w v (s) + 1 £ Peq {eM) W;( S )E 3 (0)W v (s) I . 

We will separately take the double limit as in (|1.27p for both above terms. 



2.1 The contribution of the discrete spectrum 

Let us start our analysis with the pure-point part and compute 



s — lim 

■n\o 



n-]imW v (s)*E j (0)Wr l (s) 

s\[— oo 



As V is a bounded analytic perturbation of H , the map [0, 1] 9 k i— > Ej (n) is - in particular 
Lipschitz continuous in the uniform topology. Thus there exists a constant C > such that: 



\WZ(8)Ej(0)Wr,(a) - W;(s)E j ( x (r ] s))W v ( S )\\ < C X {w 



s < 0. 



(2.2) 



Thus we can replace Ej(0) with the analytically continued projection Ej(x(v s )) an d the limit does 
not change. We will prove the following result (a weaker version of the gap-less adiabatic theorem, 
see [36] and references therein): 
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Proposition 2.7. Under our Hypothesis \1.2\ the following limit exists in the uniform topology and 
we have the equality: 



n — lira 



n-]imW;(s)E j ( X (vs))W v (s) 

s\, — oo 



which combined with 
Corollary 2.8. 



immediately gives: 



n — lim 



n-]imW*(s)E j (0)W v (s) 

s\j — OO 



(2.3) 



n — lim 



n - lim e lsH E ac (H)W v {s)E pp (K{l)) 

s~\ — CO 



= 0. 



(2.4) 



While (|2.3|) concludes the proof of the adiabatic limit for the discrete part of the spectrum 
(even in the uniform topology) , the limit in (|2.4[) is a technical result which will play a role in the 
contribution of the continuous spectrum. 



2.2 The contribution of the continuous spectrum 

We will now focus our attention on the term coming from the absolutely continuous part of the 
spectrum: 

(2.5) 



s — lim 

r,\0 



s - lim W;(s)p eq (H)E ac (H)W n (s) 

S~\ — OO 



Due to (|2.4[) we may conclude that 



s — lim 



s - lwaW*(s)p eq (H)E ac {H)W 7l {s) 

s\[— oo 



s — lim 
s — lim 



s - limW* (s)E ac {H)e- lsH Peq (H)e lsH E ac (H)W v (s) 



sNj — oo 



s - limE ac (K(l))W;(s)p eq (H)E ac (H)W r) (s)E ac (K(l)) 

s\j — CO 



(2.6) 



provided that the last double strong limit exists. Note that all errors go to zero in the uniform 
norm. 

o o 

The next step in the proof is to replace p eq (H) with p eq (H)E ac (H) in (|2.6p . In order to show 
that we can do that replacement, let us write the identity: 

{p eq {H)E &c {H) - p eq {H)}E &c {H)W n {s) (2.7) 
= -p eq (H)E pp (H)e- lsH E ac (H) {e lsH W ri {s)} + {p eq (H) - p eq (H)} e - lsH E ac (H) {e tsH W„(s)} . 

When s — > — oo both terms on the right hand side converge to zero due to the fact that e lsH W v (s) 
is convergent in the operator norm, Ce~ ltA P ac (A) converges strongly to for any A sclfadjoint and 

o o 

C compact [371 Lem.1,1 §4.4] and using the fact that p eq (H)E pp (H) is compact and the following- 
result (see § 15.11 for the proof): 

Proposition 2.9. For any continuous function $ G C(R) which tends to zero to infinity, we have 

o 

that <&(H) — $(i?) is a compact operator. 
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Up to now we have shown that the limit in (|2.6[) must equal 



s - lim i s - lim E a . c (K(l))W:(s)E a , c (H)p eq (H)E a , c (H)E ajC (H)W v (s)E a , c (K(l)) \ . (2.8) 

r;\() s\ — oo J 

For the next step we will need a comparison dynamics for W v (t), generated by the operator with 
internal Dirichlet walls. To the decoupled Hamiltonian we can associate: 

K( X (vt)) :=°H + x(vt)V. (2.9) 

o 

The associated evolution W v (t) is defined as the solution of the following Cauchy problem: 

id t W v (t) = K v (t)W n {t) 
W v (0) = 1 

(its existence results by arguments similar to those concerning the existence of W v (t)). 

o 

An important observation is the fact that A^j commutes with V so that we have 

Wr,(t) = e~ u " 1 + n_ (e~ iv - Si x{nu)du - l) + n+ (e~ iv + So x(v»)du _ {j ( 2 

with the exponentials in the second factor being just complex numbers. All terms commute which 
each other. Therefore the limit in (|2.8[) must equal: 



s - lim i s - limS ac (^(l))W';( S )S ac (F)^( S )p eg (iJ) J B ac (if)W l; (s) J B ac (F)W' ?; (s)i; ac (X(l)) 

T)\0 { S\~ OO 

(2.11) 

We state a result which will be proved later ( see § 15. ip : 
Proposition 2.10. The following limits exist in the strong operator topology: 

3,:= lim E &c {K{l))W*Js)E &c {H)W v {s)E &c {H). (2.12) 
One can see that the product of operators in the limit (|2.12l) coincides with the product of 

o o 

operators placed at the left of p eq (H) in (|2.11[) . At the same time, at the right of p eq (H) is the 
adjoint of the same product. 

Now if we can prove that S* can be written in the following way: 

E* = s - lim E ac (F)W*(s) J B ac ( J ff)W^(s)Kc(i<'(l)), (2.13) 

sNj— oo 

then the limit s — > —oo in (|2.1ip would give: 

Z v p(H)~;. (2.14) 
Indeed, since Proposition 12. 101 implies the existence of the weak limit: 

E* = w - lim E &c {H)W v {s)E &c (H)W v {s)E &c {K(l)), 

s\, — oo 

then (|2.13[) holds if we can prove the existence of a strong limit. Now in order to prove that a 
strong limit exists, let us insert some operators in the following way: 

E, c (H)W ri (s)E &c (H)W ri (s)E ac (K(l)) 

= E ac (H)W v (s)e- lsH E ac (H) {e lsH W v {s)} E ac (K(l)) 

= E ac (H) jwy S )e-^ J ^e is ° H e- isH E ac {H)} {e* sH W v (s)} E ac {K{l)). (2.15) 
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Let us investigate each curly bracket. The couple e lsH W v (s) converges in norm to o>* when s — > 
—00. The factor W rj (s)e~ lsH converges in norm too, sec (|2.10[) . Finally, the factor e lsH e~ tsH E ac (H) 

o 

converges strongly to the wave operator associated to the pair of Hamiltonians {H, H} as stated 
by the following proposition which we will prove later: 

o 

Proposition 2.11. The wave operator ui- := s — \ime tsH e~ lsH E ac (H) exists as a unitary map 

— 00 

o 

from EaciliyH. onto E ac {H)'H and one has: 

s - \ime lsH e- ls ° H E ac {H) = lo*_ = E ac (H)u>*_. 

s~\— 00 

Now we can introduce (|2.12l) and (|2.13l) in (|2.11[) , and see that the contribution coming from 
the continuous part of the spectrum will be: 

s - lim E n p(H)E* (2.16) 
The next step in our strategy is to prove that the strong limits of S^ and S* exist when 77 \ 0, 

o 

and they will equal the wave operators associated to the pair of Hamiltonians {K(l), K(l)}. First, 
we need to be sure that these limiting operators exist and are complete, and this is stated by the 
following proposition: 
Proposition 2.12. 

o o 

1. For any n £ [0, 1] we have E ac (K (n)) = E ac {H). 

2. The following limits exist: 

s - lime^We-^W^^) = : S - E ac (K(l))E E ac (H); 

s\j— OO 

s - \\me isk ^e- i8K ^E ac {K{l)) = -*= E ac (H)Z*E ac (K(l j). (2.17) 

o 

Thus the wave operators associated to the pair {K(l),K(l)} exist and are complete. 



The next technical result establishes the adiabatic limit for the wave operators S^; note that 
Dollard |16j investigated a related problem in the case of short range and relatively bounded 
perturbations. 

Proposition 2.13. S„ has a strong limit when rj \, and moreover s — lim S„ = So, where So 

O 

is the stationary wave operator associated to the pair {i^(l), K(l)} and is unitary as a map from 
E ac {H) onto E ac {K{\)).. 

We see that the very last thing to be shown in order to finish the computation of the adiabatic 
limit in (|2.16|) . is the strong convergence of S* to Sq when 7/ \ 0. Due to the completeness of the 

o 

wave operator S (point (2) in Proposition l2.12j) . we have that Sg : E ac (K(l))H —> Ea^H)^ is a 
unitary operator. Then: 




for any / £ E ac (K(l))'H and thus we have strong convergence of ^* to Sq when 77 \ on 
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With this, the proof of the adiabatic limit in (I1.19[) is concluded. 



The next sections of the paper are devoted to the proofs of the above stated Propositions 12 . 71 
[213] and Corollary [ 



3 Absence of singular continuous spectrum 

We give here the proof of the absence of the singular continuous spectrum for K(k) by establishing 
a limiting absorption principle. The main technical result of this section is the following lemma: 

Lemma 3.14. Let k <E [0, 1]. There exists a discrete set 91 C R such that for any closed interval 
I C R+ \ 01 we have the estimate (here (x) := \/x 2 + I): 



sup 

z£{x+iy\xel ,0<y<6} 



-{Q 



l} R K (z)e~ {Ql) <C(I,6,k) 



< 00. 



(3.1) 



In particular, K(k) has no singular continuous spectrum. 



Proof. We use geometric perturbation theory. Let us define a quadratic partition of unity in the 
following way: 

X- + xl + x\ = 1, X± e C°°(R), X ±(x) = 1 for ± x > 2a, X± (x) = for \x\ < a 

XoeC°°(R), X o(x) = for |a;| > 2a, X o{x) = 1 for \x\ < a. 

Fix some L > 2a. Introduce the operator K K x obtained from K(k) on the region L n (-L, L) by 
imposing Dirichlet boundary conditions at x = ±L. The operator K k l has compact resolvent, 
and let us denote it with R K x(z). Here z e C \ a(K K x). 

Now let us define an approximation for R K {z) by the following formula: 

R K (z) := x-(Qi)L(z)x-(Qi) + Xo(Qi)RkMxo(Qi) + X+(Qi)L(z)x+(Qi)- 

o 

Note that on the support of X±(Qi) the differential operators K(k) and K(k) coincide, while on 
the support of XoiQi) the operators K(k) and K K x coincide, so that we can write 



H,X-(Qi) 



R k (z)x-(Qi) + [H,x (Qi)]R k ,l( z )Xo(Qi 



R K (z) X +(Qi). 



{K(k)-z)R k (z)=U + 

+ H, x+ (Qi) 

The above commutators are first order differential operators: 

H, X ±(Qi) =-2i X ±(Qi)Pi-X±(Qi), 
[H, X o(Qi)] = -2i X ' (Qi)Pi - Xo(Qi). 
can be put in the following form: 

(K{k) - z)R K {z) = Id + X(z), e^X{z) e M{H), 



(3.2) 



Thus 



(3.3) 



(3.4) 



where the boundedness of e^ 1 ' ' X(z) is due to the compact support of the functions appearing on 
the left-hand side of the operator X(z). Thus we can write the identity: 



-(Q 



l} R K (z)e~ {Ql) = e~ {Ql) R K (z)e- {Ql) - e~ {Ql) R K (z)X{z)e- {Ql) . 



11 



Since for large values of Im (z) the norm of e^^X(z) tends to 0, we can write at least for those 
values of z that: 



5-( Q ^R K (z) e - {Ql} = e - {Ql) R K (z)e-( Q ^ \l + e {Ql) X(z)e~ {Ql) 



Now e^ 1 ^ X{z)e~^^ is compact and analytic in the upper complex plane, and has a bounded 
limit from above on any interval / which avoids the discrete set of thresholds in the leads and the 
discrete spectrum of K Rt L- Moreover, due to the exponential decaying weight on the right and 
the compactly supported cut-offs on the left, eS®^ X(z)e~^ 1 ^ can be analytically continued to 
the set {x + iy\x <G I, —5 < y < 5} for S small enough. Thus we can apply the analytic Fredholm 

alternative on this set and conclude that 
set of points. 



[l + e^ 1 ^ X(z)e (Q 1 }] exists on / outside a discrete 



4 Adiabatic limit of the discrete subspace 

In order to simplify our presentation, we adopt the conditions of Hypothesis 11.31 which means that 
we have N = 2 discrete eigenvalues which might cross at only one point when k varies. Moreover, 
they remain well isolated from the continuous spectrum. Under these conditions, Rellich's Theorem 
(Theorem 11.61 in [50)) states that the two eigenvalues are given by two real analytic functions 
{£j(K)}j g { li2 } defined for k e [0, 1]. If they cross at n € (0, 1) and only there, then there must 
exist two constants C > 0, M S N* such that 



|£i(k)-£2(k)| > C\k-k q \ m , k€ [0,1] 



(4.1) 



Moreover, their corresponding orthogonal projections Ej(k) can also be chosen to be real analytic 
on [0,1]. 



4.1 Proof of Proposition [2771 

Let us focus on j = 1. We will have to show the equality: 



n — lim 



n — lim_B ?? (s) 



with B v (s) := WVs)*£i(xM)W^(s). 



(4.2) 



This follows clearly from the next result. 
Lemma 4.15. We have: 



B v (0) =E 1 ( X {0)) =E 1 (1) and lim { sup \\B v {s) - B n 

n\o l s <o 



(4.3) 



Proof. The first two equalitites are obvious. For the limit let us remember that there exists a 
unique critical time to < when x(^o) = K o which corresponds to the intersection of the two 
eigenvalues. Fix some < 6 < 1 (to be chosen later in a more precise way). We split the negative 
semi-axis R_ in three parts: 



u 



to 



'/ 



to + rf 



,0 



(4.4) 



In what follows we will investigate how B v (-) changes when s goes through each sub-interval. 

This is the "gap- 



Near the crossing: Let us first consider the interval in the middle 

less region", but nevertheless, it is easiest to deal with. From the definition of B v (s) in (|4.2[) . and 
since K(k) commutes with E\(k), we have the important identity: 



d s B v ( s ) = V x , (vs)W;(s)E' 1 ( X (r,s))Wr l (s), 



(4.5) 
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where E[(k) is uniformly bounded in k € [0, 1] due to the real analyticity of the projector. We 
write: 



B„ 



This implies: 



to + rf 



B v 



— Br, 



to + rf 



to~V d 



— Br, 



tg±n° 



d s B v (s)ds. 



to - rf 



< Crf. 



(4.6) 



(4.7) 



Outside the crossing: In the other two intervals the eigenvalue £i(x( 7 7 s )) is isolated from the rest 
of the spectrum, as can be inferred from our Hypothesis 11.21 and 11.31 More precisely, let us show 
that it is situated at a distance larger than Crj MS than the rest of the spectrum. Indeed, using 
the splitting from (|4.1[) we may write 



|ei(xM) " ea(xM)l > C|xM - x(*o)| A/ > Cr, MS 

for every s situated at a distance larger than r/~ 1+s from to/77. Here C > can be chosen uniformly 
in 77 because we assumed that x'(t) > and \x"\ is intcgrablc. 

It means that we can find a positively oriented simple contour which only contains £i(x(v s )) 
and the following estimate holds true: 



Dr, '■= SUp 

t -n 5 to+^i] z ^ r v 



sup \\(K(x(vs)) - zy'w < c v 



-MS 



(4.8) 



sGR_\ 



We can choose the length of the contour T v to be of order 1/D V . We will treat this region by using 
a second order adiabatic development for the quasi-eigenprojector given by the adiabatic theory 
(see [29l [36] and references therein). If 



X(k) := [Ei (k)E[ (k)Ei (k) - E 1 {k)E[( K )E^{k)] 



we define: 



F v (s) := B ri (s)+r ] x , (vs)W*(s)Y(x(r)s))W ri (s) 

Y( X (vs)) ■= ~ I dz (K( X (r,s)) - z)' 1 X( X (vs))(K( X (vs)) - z)~\ 



(4.9) 



where the operator Y{n) is a solution to the commutator equation i[K(n), Y(k)] = —E[(k). The 
operator F v (s) is constructed in such way that when we compute d s F r) (s), the term 9 s B^(s) gets 
canceled and we have the identity: 



d s F v (s) 



(4.10) 



v 2 w;(s)\d x ^ I dz(K( X (x))-z) "xixix^Kixix^-z) 1 



Note that Y(x(rjs)) is a bounded operator obeying 



Wr,(a). 



X — TjS 



\Y( X (vs))\\ < Cr, 



-MS 



to-rj 6 t Q + rf 



(4.11) 



which is a consequence of 
write a rough bound of the type 



and because our choice of the contour r ?) . It follows that we can 



||0.F„OO|| < Crf- 2M \\ X '(vs)\ 2 + |x"M|), »6E-\ 



to -V S t + v 5 



(4.12) 
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Thus on any sub-interval [si, s?] of the negative real axis where the above estimate holds true we 
can write: 

||^( Sl )-i^( S2 )||<CV- 2M *, (4.13) 

due to the intcgrability properties of \ ( see ljl-10[) ). From (|4.9[) and (|4.1ip we can derive the 
estimate: 

||B,( S )-F,( S )||<Cy-^, SG E_\[^-^,^±^1. (4.14) 



V 



V 



Up to a use of the triangle inequality, on any sub-interval [si,S2] of R_ \ 
write: 



tn—y tp+r) 



\\B v { Sl ) - <Cr^ 2M . 
This estimate together with ()4.7|) imply: 

\\B v {s) - B„(0)|| = \\W:(s)E 1 ( X (vs))W rt (s) - E 1 (l)\\ < C(rf + r^ 15 ) 



we can 



(4.15) 



s < 0. 



Choose now any S <G (0, 1/(2M)); then (|4.3[) is proved, which concludes the lemma. ■ 
4.2 Proof of Corollary [2781 

The limit in (|2.3p is a trivial consequence of (|2.2[) and the result of Proposition 12. 71 The proof of 
(|2.4j) is a bit longer. We start with a lemma: 



Lemma 4.16. At fixed r\ > 0, we have the limit n — lim£? l; (s) = uj^Ei^oj*. 

s\(— oo 

Proof. This can be seen by writing 

B«(s) = {W^^e-^Ue^E^xivs^e-^Ue^W^)} ^(0)^ 

where we used the fact that each parenthesis converges in norm (even though not uniformly in rj). 



Corollary 4.17. (of Lemmas \4-15\ and \4-l(>\ ) 



lim ||w_Si(0)w*-^i(l)|| =0. 

ri\0 



(4.16) 



We can now prove (12.41) : 



n — lim 

r)\0 



n - \ime tsH E ac {H)Wr,{s)E pp (K{l)) 

s\, — oo 



= n — lim 

r,\0 

= n — lim 

»7\,0 



n - lim {e isH W v (s)} Wl{s){l - E W {H)}W V (s)E pp (K{l)) 

s\j — OO 



u>* n - lim W;(s){l - E pp (H)}W v (s)E pp (K(l)) 

s\j — OO 



n - lim [w* {E pp (K(l)) - ujr,E„(K{0)) U ;}] E pp (K(l)) = 0, 

>7\i0 



(4.17) 



where we used Corollary |4.17[ and the fact that oj* = n — \\me lsH W v {s) has norm one. 

sNj — oo 



5 Existence and completeness of stationary wave operators 



In this Section we analyze the pair of Hamiltonians {H + kV, H + nV} and prove Propositions 
EM ETUI and [231 
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5.1 Proof of Proposition 12.91 

We can approximate the function $ in the uniform norm with a sequence of Cq°(R) functions 
$ n . Thus if we can prove the proposition for smooth and compactly supported functions, then 
we are done. For such a <!>„ we can apply for example the Helffer-Sjostrand formula (or any other 
norm convergent functional calculus involving the resolvent) and argue that we can approximate 

o 

in norm the difference — $ n (_ff) with a linear combination of differences of resolvents of 

the type 

N 



Yu°i {{H-z^-iH-Zj)- 1 , (5.1) 

.7 = 1 ' 



J-- 

where Cj are complex coefficients and zj are complex numbers with nonzero imaginary part. Thus 
one can reduce the problem to showing that 

(H - z)- 1 — (H - z)- 1 =: R(z) - R(z) 
is compact for some z with Im (z) > 0. 

o 

Our decoupled Hamiltonian H (sec (|1.16|1 - ||1.18|) ) is a direct sum of three commuting operators, 

o 

and we have a sc (H) = and 

a pp (H) = a pp (U HU ) = a(U HU ) C R+. 

Let us denote by {w„}„ e N the complete orthonormal set of eigenvectors of £x> in L 2 (T>) (see 
(|1.17[l ). having eigenvalues {A n }„ S N so that <j pp (2,t>) = {A n }„gN; let Pn be the 1-dimensional 
orthogonal projection on w n in L 2 (X>). In particular, 

a ac (H) =cr ac (n_ifn_ ®n+i/n + ) = (7(n_^n_ ® n+An+) = [Ai,oo). 

Then for z E C \ [0, oo) we have 

R{z)= © [([_ - (z - Xn))"^- © (1+ - (z- X n )y 1 TT+] P n ® {-A D +w- z)" 1 

n€N L J 

with ir± : L 2 (E) — > L 2 (I±) the orthogonal projections and (l — (z — A n )) 1 the resolvent of the 
longitudinal kinetic energy on X± with Dirichlct conditions at ±o. 

o 

In order to study the Hamiltonian H and its relation with H, let us first observe that they are 
two self-adjoint extensions of the same symmetric operator 

D Q := -A D + w : C °°(L_ UCU L + ) -> H. 

Let Dq be the adjoint of this symmetric operator. In order to compare the two resolvents, R(z) 

o o 

and R{z) for z€C\[0, oo), we note that Dq extends both self-adjoint operators H and H so that: 

R(z) - R(z)) H C ker (D* - z) . 



Notice that for u £ ker (Dq — z) the distribution DqU — zu has support in the part of the boundary 
T> U X> + , where ~D± := L± n ({±a} x E d ); thus on L_ U L + they satisfy the equation: 

o 

— A.d i ±u± = zu± 

with the boundary condition u±| 2 - ±xax) = 0, for u± := . Then standard arguments show 
that our vectors u± G T-L± must be of the form u± — © a n (u± n © w n ), where the functions 

nGN 
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u±, n G L 2 (X±) satisfy the equation \±u± :Tl = (z — \ n )u±^ n . Thus w±.„ = (3±,n^ ±:nX with C±,n 
the unique complex square root of z — A„ having ±i?e£±.n < 0. Let us observe that due to the 
fact that z6C\ [0, oo) and A„ > 0, our sequence {|i?e£± ira |} contains strictly positive numbers, 
and moreover, diverges with n. Thus the infimum below is positive: 



7oO) := inf \ReC±, n \ > 0. 



(5.2) 



We have thus proved the following statement (here Qi is the multiplication operator with the 
longitudinal coordinate) : 

Lemma 5.18. Let z G C \ [0, oo) and j± G M + \ {0} be such that < j± < ~fo{z) (defined in \5.ty) , 
then we have the following estimations of exponential decay: 



e ±7±Qi n± (i?( z ) - i?(z))e ±7±Qi n± 



(5.3) 



and for ^ a ( x ) '■= e 



(with < a < 'Jq(z)) we have: 



* a {Qi)(R{z)-R(z))* a {Qi 



< c. 



(5.4) 



Taking into account that R{z)% and R{z)'H are both contained in H 1 (L^) ® H 1 (C) ® H 1 (L + ) , the 
above estimate ()5.4j) together with the compactness of Sobolev embeddings for compact domains, 



imply that R(z) — R{z) are compact operators for any z£C\ 



5.2 Proof of Proposition 12.111 

The absence of singular continuous spectrum will be proved later on in Lemma 13.14) here we only 
show completeness of wave operators by using the Birman-Kuroda method |37j . In other words, 
we want to show that the difference between some large enough powers of the resolvents is a trace 
class operator. 

We need to elaborate on the previous definition of ^> a £ C°°(R) introduced in Lemma [5.181 
We now allow any ael and further more: 

* a (») = e aV * T + I , so that \(d s \n^ a )(x)\ < C s \a\, s > 1, x € K. (5.5) 

We observe that v I / Q ,(x) is invertible everywhere on K, and if a > then VP" 1 € L fc (R) for any 
k > 1. Another fact wc shall use here is that the following multiple commutator is bounded: 

[Pi, [Pi,... [Pi, w}...]] (5.6) 

where Pi := —id x and x denotes the first (longitudinal) variable of 

Lemma 5.19. Fix z G C \ [0, oo). Then there exist kd G N large enough and ot(z) > small 

o o 

enough such that for any k > kd the operators R{z) k and ^ a {Qi)R{z) k ^! aiQi) 1 with \a\ < a(z) 
are bounded operators from L 2 (L) into BC(L) (the bounded continuous functions on L). 

o 

Proof. The result for R(z) k is based on the fact that 

R(z) k n C (lf 2fe (L_)f|^o 1 (L-)) © (# 2fc (C)fK(C)) © (# 2fc (t + )fK(t + )) (5-7) 

which is based on the commutator estimate in (|5.6j) . Now if k is large enough (depending only on 
the dimension d), the right hand side becomes a subset of BC(L). 
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In order to prove the same inclusion for the other operator, let us note that we can use a 

o o 

Combes-Thomas type rotation [TU]: ^ a (Qi)H lt i/ a (Qi)^ 1 = H + T a , where T a is a first order 
differential operator which has the following mapping property: 

T a : H k (L_) © H k (C) © H k (L+) — > if fe " 1 (L_) © if fe - 1 (C) © F fc " 1 (L + ). (5.8) 

Now if |a| is small enough, one can prove by induction with respect to k that 

k 

\\H * a {Qi)R(z) k ^ a {Qx)- l \\ <oo, (5.9) 

which means: 

* Q (g 1 ) J R(z) fc * Q (Q 1 )- 1 H C H 2k (L_) © ff 2fe (C) © H 2k (L + ) (5.10) 
and we are done. ■ 



Corollary 5.20. Let F G L 2 (R). TTien £/iere exists ^eN depending on the dimension d such that 

o o 

for any z G C\[0,oo) and any k > k d we have that F{Q 1 )R{z) k and F(Qi)^ a (Qi)R(z) k ^ a (Q 1 y 1 
with \a\ < a{z), are HUbert- Schmidt operators on %. 

o o 

Proof. Let us denote by T either R{z) k or ^ a (Qi)R(z) kl i l a (Qi)^ 1 appearing in the previous 
lemma. For any f £H we have that Tf is a bounded and continuous function. Then for any fixed 
x G L, the mapping 

n 3 f ^ (T/)(x) G C (5.11) 

defines a bounded linear functional on H, uniformly bounded in x. The Ricsz representation 
theorem allows us to conclude that T has an integral kernel obeying 

sup / |T(x,y)| 2 dy < oo. (5.12) 
xeL-'L 

Thus for any function F G L 2 (M), the operators F(Qi)T have integral kernels of class L 2 (L x L), 
hence they are Hilbert-Schmidt. ■ 



Lemma 5.21. Fix z G C \ [0, oo). Then there exist kd G N large enough and a{z) > small 
enough, such that for any k > kd we have that F(Qi)R(z) k and F(Qi) 1 $ a (Qi)R(z) k '^ a (Qi)~ 1 
with \a\ < a(z) are Hilbert-Schmidt operators on T~L for any measurable function F G L 2 (IR). 

o 

Proof. We use similar arguments with those for H but this time repeated for H. We do not give 
further details. ■ 



The final technical result needed for the Birman-Kuroda theorem is the following: 
Lemma 5.22. Let z G C\ [0, oo). Then there exists njCN large enough such that for any n > n d 

o 

we have that [R(z) n — i?(z) n ] G Bi('H) (the set of trace class operators) . 

Proof. Let us fix z G C \ [0, oo). We start with the formula (valid for any p G N): 

R(z) p - R(z) p = R{zy(R{z)-R{z))R(z) p - 1 - j . (5.13) 

0<j<p-l 

Let us choose < a < min{a(z), 70 (z)}, where 70(2) is the same as in Lemma l5.18[ Let us choose 
P > 4/^ + 1 and observe that in (|5.13[) we cither have j > 2kd or p — j — 1 > 2fc^. The idea is to 
prove that operators of the type 

Rizf^dQ!)- 1 or ^(Qi)" 1 ^) 2 *- 
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are trace class, which together with Lemma |5 . 1 81 would finish the proof. Indeed, we can write: 
fl(z) a ***a(Qi) _1 = {*_ a/a (Qi)[* a/2 (Qi)i2(«)***_ a/2 (Qi)]} 

• {^- a /2(Ql)[^ a (Ql)R(z) kd ^- a (Ql)}} , (5.14) 

where the right hand side is - according to Lemma 15.211 - a product of two Hilbcrt-Schmidt 
operators. The other operator can be treated in a similar way, up to taking the adjoint. The proof 
is over. ■ 



5.3 Proof of Proposition I27T21 

o o 

We now want to study the pair of Hamiltonians K(k) = H + kV and K(k) = H + kV, for any 
k £ [0, 1] and prove Proposition ^. 121 The only difficulty comes from the fact that the perturbation 
V has a singular commutator with H and, at the same time, it does not tend to zero at infinity. 

o o 

But for K(k) there is no difficulty due to the fact that H commutes with the bias V, while 
the last one is just a multiple of the identity on each orthogonal subspace in the decomposition 
H = H-®H ®H+. 

In fact, the only result which cannot be obtained just like in the previous section is the following 
lemma: 

Lemma 5.23. Fix z £ C \ [0, oo). There exists k' d £ N large enough such that for any any k > k' d 
we have that the operators F{Qi)R K {z) k and F(Qi)^ a (Qi)R K (z) ^ a (Qi) with\a\ < a(z) (here 
^ a is as in i5.5\) ) are Hilbert- Schmidt operators on % for any function F £ i 2 (M). 

Proof. The perturbation V is still relatively bounded with respect to H with relative bound, but 
its commutator with H defined as a scsquilinear form on the domain of H is singular. We observe 
that the main difficulty comes from the fact that the range of the operator R K (z) k is no longer 
contained in the Sobolev space H 2k (L) but only in H 2 (L) for any k £ N, due to the singularity of 
the commutator of —A with V. In fact the situation is a bit better due to the fact that V (being 
constant in the X>-space) commutes with all derivatives with respect to directions from T>. Thus, 
using the results in [24] . we may conclude that: 

R K (z) k H C H 2 (R;H 2( - k - 1) (R d )) nL 2 (L) C BC(R; H 2 ^'^ (R d )) n L 2 (h) c 

C BC(K; BC(R d )) n £ 2 (L) c BC(L) 

for k > k' d depending only on the dimension d. Thus the proof goes on exactly as in Section [5] and 
we are done with F{Qi)R K {z) k . 

Regarding F{Qi)^ a {Qi)R K {z) k ^ aiQi) -1 , we need to replace ^ a with a function Vtcwhich is 
constant in a small neighborhood of ±a. In this case, when we write 

* a K{K)*- 1 =K{n)+f a 

we see that T a equals zero around the points where V is discontinuous. Therefore, if |a| is small 
enough we will have 

WKin^Mz)^-^ < co, (5.15) 

and the proof goes in the same way as in the previous section. ■ 

6 Study of and its adiabatic limit 

Let us recall a few facts about the decoupled system: 

a sc (K( K )) = 0, H ac (K( K )) = H- ffi U+, U pp (K(k)) = U , Vac 6 [0, 1], 
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H 



U ac (H) 

n 



n_ [( - Az,._)]n_ © n+[( - A D ,+)]n+ 



w v (s)E ac (H) = [n_ ( e - i »-/o , '<w <l j + n+ U-to+JS x{vu)du\ 



-isH 



E ac (H) 



and using the notations defined earlier (|1.17[) , it is well known that: 

A D ,± = [±® l + l®£x>, cr(I±) = cr ac (l±) = [0,00), cr(£x>) = cr pp (£x.) C K+. 

o 

Thus cr ac (H) = [inf er(£x>), oo) and has the set of thresholds T = er pp (£x>)- 

6.1 Proof of Proposition 12.101 

Here we are interested in the strong limit when s — > — oo of: 

E ac {K 1 )W;(s)E^(H)W n (s)E ac (H). 

Let us start by noting that we can replace E ac (H) with the identity in the above product, and 
still get the same strong limit (if it exists). The explanation is that we can write: 

E pp (H)W n (s)E ac (H) = {E w {H)e- ls ° H E ac (H)} e isk W v {s) 

o o O O 

and use the fact that e lsH W r) (s) converges in norm, while E pp (H)e~' lsH E ac (H) converges strongly 
to zero when s — > — oo because E pp {H) is compact. Thus it is enough to study the existence of a 
strong limit when s — > — oo of: 

S„(.s) := E^K^W^sW^E^H). 

O 

For any 8 > let Vs be the set of vectors / £ % ac (H) with compact spectral support with 

o o 

respect to H at distance larger than 8 from all thresholds. Clearly, {Vs}s>{) is dense in T-L ac {K K ) = 

o 

% ac {H). It is thus enough to show the existence of the limit lim B, TJ (s)f for / £ V$. As any 

sNj — oo 

o 

vector / £ Hac(H) is of the form (/_, /+) £ %- © we will treat the two situations separately. 
The idea is to use a variant of Cook's method. We have the following identities: 

S r ,( S )/ - E ac (K 1 )W;(s)(K( X (vs)) + iy 2 W v (s)(K( x (vs)) + iff 

= E ac (K 1 )W;(s)(K( X (vs)) + l) _1 (A-(xM) + l)"^,( S )(i(xW) + iff 



E^K^W^s) 



(K( X (t]s)) + 1) 1 -(k( x ( V s)) + l) 1 



Wr,(s)(K( X (vs)) + l)f 



3.1 



Without loss of generality, let us assume that / £ H+ H Vs- Since / is with compact support in 

o 

the spectral measure of H, there exist a finite number N of transverse eigenvectors {w n } of £x> 
in L 2 (T>) corresponding to the eigenvalues {A„} and so that 



N 

/(x,x_l) = y]w n {x±) / sin[k(x - a)]f n (k)dk 
. , Jr 



3.2) 
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where f n are smooth, compactly supported, with a support which does not contain the points 
k 2 < 5. Then we have 



N 



{e- isH f}(x,x ± ) = J2™n(x±) / e 



71=1 



-is(fc 2 +A„) 



sin[fc(a; — a)]f n (k)dk, 



-is(k 2 +\ n ) 



sin[fc(x — a)]f n {k)dk. (6.3) 



{W v (s)f}(x,x ± ) = e -^Kx(nt)dt V Wn (x ± ) / 

n=i 

Moreover, for j > 1 we have: 

{I^,(s)(%W) + 1) J '/}(*,xl) (6.4) 
= e-™+ /o" ^ w;„(x_l) J e -^ 2 +*«) [fc 2 + w+x(ry.s) + A„] J sin[fc(a; - a)]f n (k)dk. 



n=l 



By standard integration by parts arguments, due to the support condition of /„, we can prove the 
following estimate: 



{W n (s)(K(x(vs)) + lYf}(x, X± ) 
or (taking N = 2) 

' -»(ft>{^W(%W) + i) J /} 



< (l + lx-oDCWC/.j), H>1,VJVGN (6.5) 



< 



C(f,j,a) 
1 + s 2 ' 



(6.6) 



We need only one more ingredient. Using the same methods as in subsection 15.11 one can prove 
the following estimation similar to (|5 .3[) for the given time dependent objects: 



e° (Ql) n±[(]f(xW) + l) 1 - (tf(xM) + 1) 'je^^IId 



(6.7) 



for a small enough. The above constant can be chosen uniformly with respect to s. 

Now we can go back to (|6.1|) and investigate the structure of ^ v (s). Let us show that it 
will converge to zero when s — > — oo. Indeed, the difference of resolvents provides the exponential 
localization near the sample. But then we know that the adiabatic decoupled free evolution decays 
with s, as in (|6.6[) . We conclude: 

Urn E v (s)f= lim $„(«), (6.8) 

s— — oo s— oo 

provided that the limit on the right hand side exists. We shall show that § v (s) has an absolutely 
integrable derivative with respect to s. Let us differentiate $ r) (s) with respect to s. We obtain 
the identity: 



id s $ v (s) = 



(6.9) 



(A'(xW) + i) 1 - (£(xM) + 1) 1 



ff,( s )(%W) + i) 2 / 



^ X '(77s)^ac(^l)^;( S )(X(xM) + 1) 



(i^(xM) + i) '-fcwi + i) 1 



W n (s){K( X (ris)) + l)f. 



We see that using (|6 . T[) and (|6.8[) we can write: 



C 



(6.10) 
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Thus lirn.s_j._oo $^(s) exists and equals: 



-i E^K^W^s) 



(j_(xM) + i) 1 -(i_(xM) + i)" 



(6.11) 

W v (s)(k( X (vs)) + l) 2 fds 



rii^E^K^Wlis^Kixiqs)) + 1) 1 



• V 



(K( X (vs)) + l) 1 -(k( X ( V s)) + l) 1 



The proof of Proposition 12.101 is over. 
6.2 Proof of Proposition 127131 



First, let us compute the limit rj \ in (|6.11|) . We can apply the Lcbcsgue dominated convergence 
theorem in (|6.11j) and obtain: 



lini iS,/ = ^(JfO^l) + 1) (K(l) + l)f 



Sac(^i)e lsK(1) 



(_:(i) + i)- 1 -(A-(i) + i)- 1 



(6.12) 

e -^M(j<:(l) + l) 2 /ds. 



Second, let us show that the above right hand side coincides with So/. Indeed, let us look at 

o 

the vector E^K^e"^ e - tsK ^ / ', where / g Vj. As before, we can decompose the vector as : 
$o(s) - *o(s) where 

_■„(*) := £ ac (^i)e^ (1) (^(l) + l) _1 (/v(l) + lJ-'e-^W (__(!) + iff, 



*o(«) :=^ac(^i)e 



isK~(l) 



(#(1) + 1) 1 -(__(!) + !)■ 



e" l;iK(1) (A"(l) + 1)/. 



Using the previous propagation estimates which were shown to be uniform in ?y, we can repeat 
the same argument which led us to (|6 . 1 1 1) but with rj = from the beginning. This will give a 
formula for So/ which will coincide with the right hand side of (|6.12j) . The proof is over. 
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